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CONTINUOUS TIME RANDOM WALKS AND THE CAUCHY 
PROBLEM FOR THE HEAT EQUATION 

HUGO AIMAR, GASTON BELTRITTI, AND IVANA GOMEZ 


Abstract. In this paper we deal with anomalous diffusions induced by Con¬ 
tinuous Time Random Walks - CTRW in R"'. A particle moves in R" in such 
a way that the probability density function u(-,t) of finding it in region Q of 
R" is given by J^u{x,t)dx. The dynamics of the diffusion is provided by a 
space time probability density J{x,t) compactly supported in {t > 0}. For t 
large enough, u must satisfy the equation u{x,t) = [(J — 5) * u] (x, t) where <5 is 
the Dirac delta in space time. We give a sense to a Cauchy type problem for a 
given initial density distribution /. We use Banach fixed point method to solve 
it, and we prove that under parabolic rescaling of J the equation tends weakly 
to the heat equation and that for particular kernels J the solutions tend to the 
corresponding temperatures when the scaling parameter approaches to zero. 


1. Introduction and statement of the results 

We shall be concerned with a probabilistic description of the motion of a particle 
in the space M". As usual we shall write to denote the set : x £ 

K" and f > 0}. Sometimes we shall also considerer the whole space time = 
{{x,t) ■. X & M" and t £ R}. The x variable is thought as a space variable, while t 
represents time. 

Let u{x, t) denote, for t fixed, the probability density of the position of the 
particle at time t. Precisely, for a given Borel set E in R." the quantity 'P{t,E) = 
Jg u(x, t)dx measures the probability of finding the particle in E at time t. 

The general problem is to find u{x, t) when the dynamics of the system is known 
and some initial state is given. 

Regarding the dynamics of the system we shall deal with anomalous diffusions. 
More precisely with continuous time random walks (CTRW). For a comprehensive 
introduction to the subject we refer to [6]. A CTRW in K” is provided by a space- 
time probability density function, the kernel, J(x, t) defined in In this model 

the particle has a probability density function u{x, t) of arrival at position x £ R" 
at time t > 0 which depends on the events of arrival at any y £ R”, sometimes only 
on the events of arrival at any y in some neighborhood of x, at any previous time 
s < t. Precisely, this dependence is given by the convolution in R”+^ of J with u 
itself. In other words, for t > 0 and x £ R" 

u{x,t) = {J * u){x,t) = // J{x — y,t — s)u{y, s)dyds. (1-1) 

JJR^+i 

The physical condition of the dependence of the current position of the particle 
only on the past (s < t) gives us the first natural condition on J, 

(Jl) suppJ CR"'*'^. 


The research was supported by CONICET, ANPGyT (MINGyT) and UNL. 

1 



2 


On the other hand, since J is a density in we must have 

(J2) J > 0, and 

(J3) J e and fjgn+i J{x,t)dxdt = 1. 

Following the notation in [6], the density function defined in R" by 

X{x) = / J{x,t)dt 

Jr 


is called the jump length probability function. Notice that from (Jl), X{x) = 
/g+ J{x,t)dt. On the other hand, the waiting time probability function is given 

by 

'''(J) = J{x, t)dx. 

Regarding the initial condition, let us first assume that the particle is localized 
at the origin of R" for t < 0. In other words u{x, t) = 6q{x) for t < 0. Hence, since 
u{x,t) for t > 0 needs to satisfy (1.1), from (.Jl) we must have that 


i{x,0)= J{x - y,-s)u{y,s)dyds 

J JR"+1 

= / ( / J{x - y,-s)u{y,s)dy) ds 

Jr- \JR’* / 

= / ( / Jix -y,-s)So{y)dy\ ds 

Jr- \Jr^ j 

= / J{x^—s)ds 

Jr- 


= \{x). 


In other words, the deterministic situation, the particle is at the origin for t < 0 
produces immediately at time t = 0 a random situation modeled precisely by the 
jump length probability function X{x) associated to the density J. 

More generally, if the position at time < < 0 of the particle distributes as indicates 
the density f{x), then u{x,0) = {X * f)ix). In this framework the basic initial 
problem we are interested in, takes the following form. Given J(x, t) and f{x), find 
u{x,t) for {x,t) € such that 


(P) 


u(x, t) = {J * u){x, t), 

fix), 


l{x, t) = 


X G R”, t > 0; 
t < 0; 
u(x, t), t>0. 

Sometimes, to emphasize the data J and / in (P)^ we shall write P( J, /) for the 
problem P and u( J, /) for its solution. 

Let us observe that the expected initial condition is attained since, taking t = 0 in 
the first equation in (P) we get u{x, 0) = {J*u){x, 0) = ff J{x — y, —s)f{y)dyds = 
(A*/)(x). 

We shall consider wide families of kernels J, but there is one, the parabolic mean 
value kernels, which plays a more significant role for our subsequent analysis. We 
shall use (for heat) to denote these special occurrences of J. Let us introduce 
the most known of these kernels (see [7] or [5]). Set W{x,t) to denote the 
Weierstrass kernel for t > 0 and x G M". Precisely yV{x,t) = (47rt)“”/^e“l^l /(4t)^ 
Set E = {(x.t) G : >V(a;,t) > 1} and 'H(a;,t) = jXEix,t)fl^. 
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As it is easy to check % satisfies properties (Jl), (J2) and (J3) stated above. 
Moreover, % satisfies also the following two properties 
(J4) has compact support in 
(J5) it is radial as a function of a: S K." for each t. 

The outstanding fact regarding T-i is given by its role in the mean value formula 
for temperatures. If v{x, t) is a solution of the heat equation ^ = Av in a domain 
n in then, for (x, t) € il and r small enough we have that v{x, t) = fj Hrix — 

y^t — s)v{y, s)dyds, where T-Lr denotes the parabolic r-mollifier of H. Precisely 

'Hrix,t) = ^„_|_2 ^ , 

with E{r) = {{x,t) G : W{x,t) > r~^}. The following figure depicts the 

support E{r) of Hr- 



Figure 1. Sets E{r) for n = 1 and r = i. 


In the sequel, for any kernel J{x,t) and r > 0 we shall write Jr{x,t) to de¬ 
note the parabolic approximation to the identity given by Jr{x,t) — ipkpzJ (f, ;^)- 
Moreover, the notation Vr{x,t) or even fr{x) for functions depending on space 
time or, only on the space variable will have always the same meaning. Precisely, 
Vr{x,t) = and/^(a;) = f x) . 

The results of this paper are in the spirit of those in [4] and [3]. Instead of 
dealing with generalization of boundary conditions, we are concerned with diffusion 
problems in the whole space K" and the initial condition is generalized. 

Let us state the main results of this paper. The first one is the weak convergence 
to the heat equation. 

Theorem 1. Assume that J{x,t) satisfies (J2), (J3), (J4) and (J5). Then, for 
each ip in the Schwartz class o/K”+^, we have 

1 

lim — - 5) * p = -I- vAp, 

i-->o at 

uniformly on where p-= — ff tJ{x,t)dxdt and u = ^// \x\^ J{x,t)dxdt. 

The second result concerns the existence of solutions for problem (P). For a given 
Lipschitz function of order 7 , / G we denote by [f]^ the corresponding 

seminorm of /. In the next statement C denotes the space of continuous functions. 
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Theorem 2. Assume that J{x,t) satisfies (Jl), (J2), (J3) and (J4)- Set a = 
sup{/3 : J{y, s)dyds < 1}. Let f € L°“(]R”) be given. Then there exists 

one and only one solution u{x,t) of (P) in the spaee (C fl If f G 

(L^ n L°°)(MA), then u{x,t)dx = f{x)dx for every t > 0. In partieular, 
if f is a density function, so is u{-,t) for every t > 0. Moreover, if f belongs to 
((-0,7 n we have that 

\u{x,t) - f{x)\ <C[f]j (1.2) 

for (x,t) G R” X [0,a] and some C which does not depend on f. 

The next result which is interesting by itself contains a maximum principle which 
shall be used in the proof of Theorem 4. Precisely, the supremum of the probability 
density function in the future of a = sup{/3 : Jiy, s)dyds < 1} coincides with 

its supremum in R" x [ 0 ,a]. 

Theorem 3. Let J be a kernel satisfying (Jl), (J2), (J3), and (J)). Let w{x,t) 
be a bounded function defined in R"'*'^ such that 

w{x,f) = JJ J{x — y,t — s)w{y,s)dyds (1-3) 

for {x,t) G R" X [a,+oo). Then, 

sup |i(;(x,t)|= sup \w{x,t)\. 

(x,i)eR"+^ (x,t)eR"x[0,a] 

Let us proceed to state the fourth result of the paper. 

Theorem 4. For each H G Jif there exists C > 0 such that for every r > 0 and 
every f G D L°°){M^) 

/) — u||^oo(R^+l) ^ C[f].yr'^ , 

where u is the temperature in given by u{x,t) = (>V(-,t) * f){x). 

Let us finally remark that in [2] the authors prove the Holder regularity for 
solutions of the master equation associated to CTRW’s. 

In Section 2 we prove the weak convergence of parabolic rescalings to a heat 
equation. In Section 3 we show existence of solution for the Cauchy nonlocal prob¬ 
lem. Section 4 is devoted to prove the maximum principle contained in Theorem 3. 
Finally, Section 5, deals with convergence of solutions of rescalings of (P) to tem¬ 
peratures. 


2. Some space time nonlocal parabolic operators and their weak 
LIMIT. Proof of Theorem 1 

For 0 < r < 1, since // J{y, s)dyds = 1, applying Taylor’s formula we get 
JJ Srix -y,t- s)(p{y, s)dyds - (p{x, t) 

= JJ Jr{x -y,t- s){(p{y, s) - ip{x, t))dyds 


Jr{x -y,t- s) 


XI t){s - f) 


dx. 


dt 
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+ 2 (y “ 2 ;, s - t)D^(p{x, t){y -x,s-tf + R{y -x,s-t) 


dyds, 


where denotes the Hessian matrix of the second derivatives of (p with respect 
to X and t and |i?(x,<)| = 0{\x\^ + t^)i . 

The last integral in the above identities can be written as the sums of the fol¬ 
lowing seven terms, 


^ ^ (yJJ ~ Xi)Jr{x -y,t- s)dyds^ , 

II = ^{x,t) f f f {s - t)Jr{x -y,t- s)dydsj , 


E 


d'^ip 


.. dxidxj 


ix,t){ - / / {yi - Xi){yj - Xj)Jr{x -y,t- s)dyds 


/F = ^ Q JJ(yi -XifJr(x -y,t- s)dyds^ 




2=1 


(ui - Xi)(s - t)Jr(x -y,t- s)dyda 


VI = -^(x,t) ( - 11 (s-t) Jr(x -y,t- s)dyds 


and 


VII = // Jr{x — y,t — s)R{y — X, s — t)dyds. 


Since for t fixed J is radial as a function of x, then /, III and V vanish. For the 
other four integrals we perform the parabolic change of variables (z, () = 
to obtain 

II =^(x,t)r^ (^- JJ C,J(z,C,)dzdC, 

= {IJJ zfJ(zX)dzdJ 


VI = 


d‘^(p 

W 


Mr' {I II 


ej(zx)dzdc 


VII = JJ J(z,C)Rirz,r'^OdzdC. 

Finally, since, as a function of r close to zero, VI and VII are of order at least r^, 
we see that 


!—>0 


!—>0 \ r‘ 


II IV 


l 2 -S)*ip] (x, t) = lini ( 7 : 5 - + — ) = p-^ix,t) + i/A(/?(x, t) 


dp 


where y and v are defined as in the statement of Theorem 1 . That thus convergence 
is uniform in follows from the fact that tp is a Schwartz function and so VI 

and VII converge to zero uniformly. 
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Lemma 5. For J — 'H we have that = —v and the limit equation in Theorem 1 
is the heat equation multiplied by a constant. 


Proof. All we need to show is that 


U{y,s)s dyds = — j j 'H{y,s)\y\‘‘ dyds. 


( 2 . 1 ) 


Let us compute both of them in terms of the Euler gamma function and the area 
surface of the unit ball of M", 5'"“^. On one hand we have that 

.|2 


JJ 'H{y:S)s dyds = 2 // ^E(i){-y,-s)^-^s dyds 


E{1) 


-dyds 


Iff . he,,,. 

4 J-^ Js(o,(2nsln(47r(-s)))i) 


f*{2ns ln(47r( —s))) 2 




A -Sho 


dadpds 


r(S'" 


rO 


4(n + 2) 

1 / n 


1 n +2 

— (2ns ln(47r(—s))) ^ ds 
J- —s 


4(n + 2) Jq t 
(T(S'"“^)n^ 


4(n + 2)2 2 7r"J^ Jo 


(^4-ln(t))) ' dt 

e-^{--^)e^de 


a{S'^ ^)n” 2 ^ 2 2 2 "+2 

4(n + 2)2^^^ {n + 2) (n + 2)^ Jo ^ ^ 


r(S'"-l) 


n 2 


...N " + 6 n + 2 

2 (n+2) 2 TT 2 


n + 4 


On the other, 

^// d2/ds= ^ JJ XE(i){-y,-s)^-^\yf dyds 

= i:ff Kdyds 


8n 

8n 


1 


£;(!) 


f [ . ^ 

JB^0,(2nsln(4ii-(-s)))5) ^ 


■dyds 


,0 


n{2ns ln(47r( —s))) 




8 n 
a(S'"-i) 


s^ 

iTT 

n—1\ 


/ 5 n-l 

n+4 


dadpds 


8n{n + 4) 7-i 
^(^"-i) n^47r 

8(n + 4)2^7r^ Jo 


(2ns ln(47r(—s))) ^ ds 


4 (t(-ln(t))) ^ dt 
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a(6’"-i)n^47r 
8 (n + 4)2^ TT^ 






cr(5'"-i)n^47r 


8(n + 4)2”2^7r 2 (^ + 2)(n + 2)"2 


f 


e-H^dC 


r(S'’^-i) 


n 2 


(n + 4)(n + 2) 2 tt 2 


n + 6 


Now, since r (0 + 1) = zV{z), we have that ;;;^r and the proof 

is complete. □ 


3. Existence of solutions for (P). Proof of Theorem 2 

Let J{x, t) be a kernel defined in space time ]R”+^ satisfying (J1), (J2), (J3) and 
(J4)- Let / G L°“(IR”) be given. Following the ideas in [4], [3] and [1] we shall solve 
(P) by iterated application of the Banach fixed point theorem. From (J3) and (J4), 
a = sup{,5 : J{x,s)dxds < 1} is positive and finite. For the first step in the 

use of the fixed point theorem in the Banach space = (C D L°°){W^ x [0, ^J) 
with the L°° norm. 

As in the statement oi (P) set 


* < 0 

1 te[o,f] 


where v G 3§i. Since v is bounded on R" x (— 00 , j] and J G the integral 




J Jw 


/ Q: 1 

x(-oo,^] 


J{x-y,t 


s)v{y,s)dyds 


is absolutely convergent for (x, f) G M" x [0, ^]. Let us prove that, as a function of 
{x,t) G R." X [0, y] the function g belongs to In fact, from the definition of g 
we see that 


\g{x,t)\ < (^jj Jdyds^ ||u||,^ < sup{||/||,^ , ||u||^}. 


Let us check the continuity of g. For h G M" and fc G R. such that {x + h,t + k) G 
(— 00 , y], we have that 


\g{x + h,t + k) — g{x,t)\ < JJ \J{x + h — y,t + k — s) — J{x — y,t — s)\\v{y, s)\dyds 

< wi(Y^|/i|^ + fc2) ||u||^, 


where uji is the modulus of continuity in of J. Hence for v G 3§i we also have 

that g G 3^1 when restricted to the strip R” x [0, ^]. 

Define Ti : —>■ by Tiv = g. 

Let us now prove that Ti is a contractive mapping in ^i. Let v and w be two 
functions in 3Si. Let {x,t) G R" x [0, ^]. Then, with 


w{x, t) = 


f{x), t < 0; 
w{x,t), tG[0,f], 
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we have that 


Tiv{x,t) - Tiw{x,t) = ^ J{x - y,t - s){v{y,s) - w{y,s))dyds 

= Jix - y,t - s){v{y,s) - w{y,s))dyds. 

J Jo<s<% 


Hence 


||Titi - Tiwll^ < ( sup // J{x-y,t-s)dyds\\\v-w\\^ 
^(a:,t)eR’»x[0,'^] •^ 0 <s <'2 

Now, from the definition of a, and (Jl) 

J{x — y,t — s)dyds = 


< 


J{z,a)dzda 
/ / J{z, a)dzda 

J J0<(T<t 


J{z,a)dzda =: r < 1. 


So that \\Tiv — < t\\v — w||oo- Hence Ti is a contractive mapping in 

So that there exists a fixed point t6i € for Ti; TiUi = ui. In other words 


ui{x,t) = // Jix -y,t- s)ui{y,s)dyds 


(3.1) 


for X G M” and 0 < t < f. 
Let us check that 


/ ui{x,t)dx = / f{x)dx 

JR" JR" 


for every 0 < t < ^, when / G L^(]R”). Since ui can be realized as the limit 
of the sequence of iterations of Ti applied to any function v G we may take 
v{x,t) = f{x) as the starting point. In doing so we see that the integral in variable 
X of \T^f{x,t)\ is less and equal to / |/| dx. In fact, from (JS), we see 

J \Tif{x,t)\dx = J jj J{x -y,t - s)f{y)dyds dx 

< J (^JJ Jix - y,t - s)dxds^ \fiy)\dy 

= J \f\dy- 

Hence, inductively, assuming J |T™/(x,t)| dx < J \ f\ dx we have 

J |T™+V(x,t)|dx = J \T,{T^f)ix,t)\dx 

= I JJ Jix-y,t- s)T:i^f{y, s)dyds 
= / JJ Jiy,t-s)T[^f{x-y,s)dyds 


dx 

dx 
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/ // J{y,t - s)f{x - y)dyds + 

J J J s<0 

+ [[ J{y,t - s)T^f{x-y,s)dyds 
J Js>0 


dx 


< 


J{y,t- s) 


s<0 


f{x - y)dx 


dyds+ 


+ Jiy,t- s) 


< J \f\dx. 


s>0 


J T^fix-y,s)dx 


dyds 


With the same arguments we can conclude that J f{x,t)dx = J f{x)dx for 
0 < t < §• The result then follows since, for f G we have that T^f tends 

to ui also in C([0, In fact, if we prove that 

ll|r™+V-rr/|||<r-|||ri/-/lll (3.2) 

where |||u||| = sup^^jg qj lk(■^^)llLl(Rn)) then T™/ is also a Cauchy sequence in 

C([0, ^],L^(IR")). Since T™/ converges uniformly to Ui we get the desired preser¬ 
vation of the integral. Let us prove (3.2). 

Let us first check that is continuous as a function of t G [0, with values 
in L^(]R") for each m. Take t and t + h two points in [0, ^]. Then 


\T^fix,t)-T^f{x,t + h)\dx 


J{x — y,t + h — s)T^~^f{y,s)dyds dx 
j j (J(z, t — s) — J{z,t + h — s))T^~^f{x — z, s)dzds dx 

< J J \ J{z,t — s) — J{z,t + h — s)\ f{x — z, s) dx^ dzds 

< J \f{x)\dx J J \ J{z,t — s) — J{z,t + h — s)\dzds, 


which tends to zero when h ^ 0 because J G L^(IR"+^). 

Similar calculations show that T^f is a Cauchy sequence in the ||| • |||. In fact, 
for t G [0, we have 


J \Tr+^f{x,t)-T^fix,t)\dx 

= I JJ J{x-y,t-s) (r^fiy, s) - T™" V(y. s)) dyds dx 

= f f f J{x-y,t-s){T^f{y,s)-T^~^f{y,s))dyds dx 
J J J0<s<t 

J{z, t — s) ^y* \T^f{x — z,s) — T^~^f{x — z,s)\ dx^ dzds 


0<s<i 
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[ [ — f |r™/(x, s) — ^ f{x, s)\dx\dzds 

j Jo<s<t \J ) 


< 


J0<s<t 

! 


sup j 

^|rr/(x,s 

^s€lO,q]J 



rr^fWl 


iii^r+v 


/0<s<^ 


hence 

<r|||r-/-Tf 

By iteration we obtain (3.2). 

Let us observe that since ui {x, t) can be obtained as the iteration of Ti starting at 
any function v in we can in particular take v as the constant function 
where s(/) = sup/ and i{f) = inf/. Then v = + fX[t<o}, so that 

*(/) <v< s{f) everywhere. From (J2) and (J3) we also have i{f) < Tiv < s{f) 
on M" X [0, ^]. The same argument shows that for every iteration of Tiv we 
have i(/) < T^v < s{f). Since ui is the uniform limit of T^v we get 

*(/) < ui{x,t) < s{f) 

on the strip R" x [0, ^]. So far we have existence and mass preservation for t G [0, f ]. 

Now proceed inductively by covering 1R+ with intervals of the type — 

The first step, i = 1 is precisely the one described above. Assume that Ui € = 

(C n L°°)(]R" X [(i — 1)^, *^]) for j = 1,...,/ have been built in such a way that 


with 


i{x,t) = jj J{x - y,t-s)ui{y,s)dyds, 
u7rT(a;,t), 




Ui{x,t), (i —l)f <t<*f. 

Moreover, Ui{x,t)dx = /j^„ f{x)dx for [i — 1)^ <t< 

i{f) <u^{x,t) < s{f) (3.3) 

for every {x,t) G R” x [(i — and Ui{x,{i — 1)^) = Ui-i{x,{i — 1)^) for 

every x. 

Define ^j+i as the space (CnL°°)(R" x [/^, (/ +1)^]) with the complete metric 
induced by the L°° norm. For v G ^j+i, define 

Tj+iv{x, t) = JJ J{x-y,t- s)v{y, s)dyds 
v(x, t) = 


with 


Uj{x,t), t<jf, 

v{x,t), /f<t<(/ + l)f. 

As in the case of i = 1, it easy to check that with {x,t) G R" x [j^, (/ + l)f], 
Tj+iv G ^j+i- Hence —>■ ^j+i- It is also easy to prove that Tj+i is 

contractive on with the same rate of contraction r obtained when i = 1. 

Also, with the same argument as in the case i = 1, with f uj(x, t)dx = / f(x)dx 
when t < we have that for t G [jf, (/ + l)f] 


/ Uj+i{x,t)dx = / f{x)dx. 

/r" Jr 
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In order to check that Uj+i(x,j^) = Uj(x,j^) we have to observe that for 
< t < {j + 1)^, the fixed point equation is given by 


Uj + i{x,t) 


J{x-y,t- s)uj+i{y, s)dyds. 


For t = j^, property (Jl) shows that the above integral only involves values of s 
which are bounded above by Then, for those values of s, u_,+i(y, s) = Uj{y, s). 
So that 

= jj J{x- y,j^ - s)uj{y,s)dyds = Uj{x,j§), 

as desired. 

Property (3.3) for i = j + 1 can be proved following the same argument used in 
the case i = 1 . Let us notice that the function u{x,t) defined in by u{x,t) = 

with j{t) the only positive integer for which {j{t) — 1 )^ < t < is 

continuous and bounded. Moreover, i{f) < u{x,t) < s{f ) for every {x,t) g 
The above remarks prove that m S .^ = (C fl and solves (P). 

In order to prove the uniqueness of the solution u let us argue as follows. Assume 
that u and u are two solutions. Then their restrictions on the strip K" x [0, 
coincide. Since the fixed point of Ti is unique and being a solution of (P) in 
R" X [0, is equivalent to be a fixed point for Ti, we see that u = u on M" x [0, ^]. 
For the next time interval [^, a] the restriction of both, u and u to this interval are 
fixed points of the same operator T 2 . Again the uniqueness for the Banach fixed 
point guarantees u = u on R” x [f, a]. Proceeding inductively we get that u = u 
everywhere. 


Let us finally prove the estimate (1.2). First we shall show that (1.2) holds 
true when {x,t) g M" x [0, ^]. This can be accomplished because the function u 
in the first time interval [0, j] coincides with ui which is provided by the Banach 
fixed point theorem and the rate of convergence can be estimated by the contraction 
constant r. We already know that t = fJ[y^ s)dyds < 1. Set to denote the 

TO-th iteration of Ti applied to the initial guess = /, then since 11 11 00 — 

t"* IIm} — u?|| , we see that 

M -L J- M CX3 ’ 

for every m = 1 , 2 ,... 

Let us now show that for (x, t) g R” x [0, there exists a constant C depending 
only on J such that ||u} — /||^ < C[f]^. In fact, 

\u\{x,t) - f{x)\ = \{Tif){x,t) - fix)\ 

= JJ 'J{x - y,t - s){f{y) - f{x))dyds 

<[fh (^JJ J{x-y,t-s)\x-yl"^ dydsj 


Hence for every m = 1, 2,... we have 

ll<-/llLoo(M»x[0.f])<^[/]7> 
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where C depends only on J. The same is true for the uniform limit ui of the 
sequence rt™, in other words 

11^1 — /llLoo^Rnx[0,^]) — ^[fh- (3-4) 

Let us now check how to get the same type of estimate for the time interval ci]- 
From the construction of u we have that on R" x [§, a], u = U 2 with 


U2{x,t) 


J{x — y,t — s)u2{y, s)dyds 


[ f{x), t < 0; 

u^{x,t) = ^ Ui{x,t), tG[0,f]; 

[ U2{x,t), t G [f ,a]. 

On R" X a] the solution U 2 is the only fixed point for the operator T 2 and, since 
the limit U 2 of iterations u™ of T 2 U 2 = u\ is independent of the starting point it®, 
let us take again w® = /. Hence ||it 2 — f\\^ < ||it 2 — /|| . Notice that if we 

write 7{y,s) = f{y)Xs<o{s) + ui{y,s)X^^^ a^{s) + f{y)X^a^^^{s) 



J{x-y,t- s)f{y,s)dyds. 


Let us finally check that the desired estimate holds for ||it 2 “ /|| ^ [fjCi]- 

Take {x,t) G R” x [§,a], then 


< 


\ul{x,t) -/(x)| = |(T 2 /)(x,t) -/(x)| 

= JJ Jix - y,t - s)f{y,s)dyds - f{x) 

[[ J{x - y,t - s)\f{y) - f{x)\dyds 
J Js<0 

^ J{x - y,t - s)\ui{y, s) - f{y)\ dyds 
'<2 

+ [[ Jix - y,t - s)\f (y) - f{x)\dyds. 

J J S<.Oi 


The first and the third terms in the right hand side of the above inequality are 
bounded by the product of the Lipy seminorm of / and a constant depending only 
on J. For the second term we use (3.4) and we are done. 


4. Maximum Principle: Proof of Theorem 3 

Recall that a = sup{/3 : J{y,s)dyds < 1}. Since the function I{f3) = 

ffs </3 JJyJs is increasing and continuous as a function of /3, a is also the infimum of 
those values of /3 for which /(/3) = 1. Moreover, from definition ofa, 0</(^)<l. 

Let tk = a + ik — 1)^, = R” x [0, tk], Sk = sup^^, \w\ for k = 1,2, _ Let us 

see that Sk = Sk-i- Let {x,t) G R” x [tk-i,tk], hence 
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S 


ik 

ik-1 

tk-2 




a 



\w{x,t)\ = 


jj j{x-y,t-s)w{y,s)dyds 


JJ J{x - y,t - s)w{y,s)dyds + JJ J{x - y,t - s)w{y,s)dydi 


tk-i<s<t 






< Sk // J{x - y,t - s)dyds + Sk -1 


J{x — y,t — s)dyds 


tk-l<S<t 


t — CX<S<tk — l 



( \ 


( \ 
r r 1 

Sk 

1 

1 

1 

+ Sk-i 

/ / J{x-y,t- s)dyds 


y t—ct<s<tk-i y 


\t—a<B<tk-i j 


= Sk - [Sk - Sk-i) \ I J{x -y,t- s)dyds 


Hence 


since 




J{x -y,t- s)dyds {Sk - Sk-i) < Sk - \w{x,t)\, 




J{x — y,t — s)dyds = // J{x — y, si)dydsi 

t — a<s<tk — i J Jt — tk — \'^s\<OL 

> [[ J{x-y,si)dydsi = l - I{^). 

J J ■^<Si<Q 

Then 

(l-/(f))(^fc-5fe_i) <5fc- sup \w\ (4.1) 

Of course Sk > Sk-i, if Sk > Sk-i, then Sk = |w| and the right hand of 

(4.1) equals zero. Since 1 — I{^) is positive we have that Sk < Sk-i, a contradiction. 


5. Convergence of solutions. Proof of Theorem 4 


Along this section we shall assume that J in problem (P) is a fixed parabolic 
rescaling of a mean value kernel H G More precisely for r > 0 we consider the 
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CTRW, with past density distribution f{x), generated by the space time probability 
density Hr{x,t) — ^), with H{x,t) satisfying (Jl) to (J5) and the mean 

value formula. The space probability density function for the localization of the 
moving particle at time t is the solution u{Hr, /)(•, t) of P{Hr, /). 

The result contained in Theorem 4 will be a consequence of the two following 
lemmas and the maximum principle contained in Theorem 3. 


Lemma 6. Let J be a kernel satisfying (Jl), (J2), (J3) and (Jf). Set a = sup{/3 : 
//s </3 < !}• S n O < 7 < l. Then there exists a 

constant C > 0 such that for every r > 0 


iiJrJ)- fWi 


^ X [OjCur^]) 


< c[f]^, 


The next lemma contains a well known result on the rate of convergence of 
temperatures to the initial condition when it belongs to a Lipschitz class. 


Lemma 7. Let f £ riL° 
solution of 




andu{x,t) = {Ant) 2 4t f{y)dy the 


— — Av 

dt ~ 

u{x,0) = f{x), 
Then there exists a constant C > 0 such that 


X £ 


\u{x,t) - f{x)\ < C[f]^t2 


for every {x, t) £ M"’ 


n+1 


Proof of Theorem 4- From Lemma 6 and Lemma 7 we have that 
sup \u{HrJ){x,t)-u{x,t)\ <C[f]^r'^. 

(£c,t)€K^ X [0,ar2] 

Now, since a = sup{,5 : H{y,s)dyds < 1} we also have that ar^ = sup{/3 : 

Hr{y,s)dyds < 1}. Hence for {x,t) e M" x (ar^,+oo) we have that the 
support of Hr{x — y,t — s) as a function of {y,s) is contained in So that 

for a temperature u defined on and t > ar^, since H £ , the mean value 

formula holds and 


u{x, t) 



y,t 


s)u{y, s)dyds 


for X £ M" and t > ar^. 

On the other hand, we also have that u{Hr, f) = Hj. * u{Hr, f) for x £ K" and 
t > ar"^, because u{Hr,f) solves P{Hr,f)- Hence, applying Theorem 3 with Hr 
instead of J, ar^ instead of a, u{Hr, f) — u instead of lu, we get 

snp \u{Hr,f){x,t) - u{x,t)\ < sup \u{Hr, f)ix,t) - u{x,t)\ < C[f]^r^, 

R”+i R’*x[0,ar2] 

as desired. □ 


Proof of Lemma 6. The result follows from (1.2) in Theorem 2 after parabolic 
rescaling. In fact, set u{J,g) to denote the solution in K”+^ of P(J,g). With 
this notation we have that 


u{Jr,f) 


u 



r 
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for each r > 0. Hence, for a: G K" and 0 < < a, from (1.2) with /(r-) instead /, 

we have 


\u{JrJ){x,t) - f{x)\ 


u 



{x,t) - f{x) 




<C[/(r.)]^ = Cr^[/]^. 


□ 


Proof of Lemma 1. For {x, t) G K."^ , since the Weierstrass kernel has integral 
equal to one, we have 

1 


\u{x,t) - f{x)\ = 


< 


< 


(47rt)"/^ , 

[fh 

(47rt)"/2 

[/I- 


e f{y)dy-f{x) 


e \f{y) - f{x)\dy 


la: — y\^ 

e |y _ x\^ dy 


' ./lU^ 


□ 
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